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ABSTRACT 

Context. The tumbling pattern of a bar is the main parameter characterising its dynamics. From numerical simulations, its evolution 
since bar formation is tightly linked to the dark halo in which the bar is formed through dynamical friction and angular momentum 
exchange. Observational measurements of the bar pattern speed with redshift can restrict models of galaxy formation and bar evolu- 
tion. 

Aims. We aim to determine, for the first time, the bar pattern speed evolution with redshift based on morphological measurements. 
Methods. We have selected a sample of 44 low inclination ringed galaxies from the SDSS and COSMOS surveys covering the redshift 
range 0< z <0.8 to investigate the evolution of the bar pattern speed. We have derived morphological ratios between the deprojected 
outer ring radius (Rring) and the bar size (Rbar)- This quantity is related to the parameter = RcR/Rbar used for classifiying bars in 
slow and fast rotators, and allow us to investigate possible differences with redshift. 

Results. We obtain a similar distribution of R at all redshifts. We do not find any systematic effect that could be forcing this result. 
Conclusions. The results obtained here are compatible with both, the bulk of the bar population (~ 70%) being fast-rotators and no 
evolution of the pattern speed with redshift. We argue that if bars are long-lasting structures, the results presented here imply that 
there has not been a substantial angular momentum exchange between the bar and halo, as predicted by numerical simulations. In 
consequence, this might imply that the discs of these high surface-brightness galaxies are maximal. 

Key words. Galaxies: high-redshift - Galaxies: evolution- Galaxies: structure- Galaxies: kinematics and dynamics 



1. Introduction 

Stellar bars are thought to be the main internal mechanism 
driving the dynamical and secular evolution of disc galaxies. 
They are able to modify the central structure and morphol- 
ogy of galaxies, mostly because they are responsible for a 
substantial redistrib ution of ma ss and angular moment um in 
the disc s (e.R., Sellwool Il98ll: IS ellwood & Wilkin sonL 11993 ; 
Pfennig er & Frie dli. 1991: lAthanas soula, 2003; Deba ttista et all 
|2006). In the last decade, numerical simulations have ad- 
dressed the importance of the transference of angular mo- 
mentum between bary onic a nd dark matter components (e.g., 
IPebattista & SellwoodlT998ll2000 ). The amount of angular mo- 
mentum exchanged is related to the specific properties of the 
galaxies, such as the bar mass, halo densit y, and halo velocity 
dispersion I Pebattista & Sellwood, 1998; AthanassoulaL 120031: 
ISellwood"& Debattista, 2006) and it takes place mainly at the 
disc r esonances ( Athanassoulai l2003l: iMartinez-Valpuesta et all 
l2006h . Recent works have shown that gas fraction can also play 
an important role in the bar-halo interplay and therefore in the 
bar e volution (Bournaud & Combes, 2002; Romano-Diaz et all 
[20091: IVilla- Vargas et al.1. 12010). Moreover, bars are efficient 
at funneling material toward the galaxy centre and possi- 
bly they influence the build ing of the stellar bulge (e.g., 
'Kormendv & Kennicutt, 2004) and the central black hole (e.g. , 
Shlosman et al, 1989; Regan & Mulchaey, 1999; Corsini et al.l . 
120031) . Peanut/boxy bulges in galaxies are also thought to 
be associated with bending instabilities and bar vertical 
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reson a nces (iBureau & Freemanl 1 19991: iDebattista et all 
'2006'; "Athanassoula. '20051: iMartinez-Valpuesta et all 
Mendez-Abreu et al, 2008). 

As already mentioned, the bar formation and destruction 
rate, as well as the morphological and dynamical changes 
suffered by the discs during their lifetime are strongly af- 
fected by the angular momentum exchange. Therefore, the 
cosmological evolution of the bar fraction can also depend 
on this effect. Observations show that bars in low redshift 
galaxies are ubiquitous, with a fraction o f -45% at visual 
wavelengths ( e.g., Marinova & Jogee, 20 071: iReese et all [20071; 
iBarazza et aE 2008; Aguerri et al., 2009; Masters et al., 201 1') 
and rising to -70% i n the near-infrared ( Knapen et al., 200Q: 
lEskridge et all l2000l: iMenendez-Delmestre et all l2007h . The 
bar fraction depends on morphological type, being lower in 
lenticular galaxies than in spirals (Marinova & Jogee 2007; 
Aguerri et al. 2009; Nair & Abraham 2010; Barway et al. 
2011; but see also Masters et al. 2011). Some recent results 
show that bar fraction is a strong function of galaxy mass 
(^Mendez-Abreu et all 120101: lNair& Abraham', '20 1 0) and color 
(Hoyle et al., 201 ij). In co ntrast, bar f raction is only barely af- 
fected by the environment (lAguerri et all l2009l: iLi et all l2009l: 
iMendez-Abreu et all[2QTQl) . 

The e volution of the bar frac tion with redshift is still a matter 
of debate. [Abraham et"an (Il999l) found that the fraction of barred 
galaxies at z > 0.5 is lower than the local fraction. However, 
other authors claim that this may be the consequence of selec- 
tion effects, due to the high angular resolution needed to find 
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bars (Elmegreen et al. 2004; but see van den Bergh 2002). To 
deal with the angular resolution problem, several studies have 
carried out this analysis using the Advanced Camera for Surveys 
(ACS) on the Hubble Space Telescope (HST). Thus, Elmegreen 
et al. (2004) and Jogee et al. (2004) found the same bar fraction 
(-40%) at redshift z = 1.1 as in the local Universe, suggesting 
that the bar dissolution cannot be common during a Hubble time 
unless the bar formation rate is comparable to the bar destruction 
rate. On the contrary, Sheth et al. (2008), in a recent study using 
images from t he Co smological Evolution Survey (COSMOS; 
IScoville et al.L l2007l) and using a larger sample than previous 
studies, found that the bar fraction at z = 0.84 is one-third of the 
present-day value. They also found a much stronger evolution 
for low mass galaxies and late-type morphological types. Part 
of the differences may be due to the selection effects and other 
systematic effects that still need to be investigated further. 

In any case, these results show that bars have been com- 
mon structural components of the discs of galaxies during the 
last 8 Gyrs. The study of their origin and evolution could be 
crucial for understanding the galaxy evolution since z = I. 
This study can be done by analysing three parameters that char- 
acterise the bars: length, strength and pattern speed. Several 
methods and techniques have been proposed in order to mea- 
sure these bar parameters. The bar length has been obtained di- 
rectly by visua l inspection on galaxy images ( Kormendv, 197' 



Martini Il995l: fM endez-Abreu et al. 2010; Nair & Abraham 
20101: iMasters et al.L 1201 ih . searching for the maximum ellip- 
ticity of the galaxy isophotes (IWozniak et al.lll995l:lLaine et al.L 
120021: iMarinova & Jogeel l2007l: lAguerri et all. l2009l). l ocating 
variations of the isophotal position angle (Sheth et al., 2003; 



Erwinl [2QQ5b . analyz i ng Fo urier moments (O uillen et a l., 1994 ; 



Aguerri et al.L l2000al l2003l) . or by photome tric decomposition 



of the surface-brightness profiles of galaxies ("Prieto et al.Vl997', 
[Aguerri et al, 2005; Laurikainen et al, 2005; Gadotti, 
20081: IWeinzirl et al.1. 120091: iGadottii 1201 ID . The resulting stud- 
ies reported that the typical bar length is about 3-4 kpc, and 



strongly correlates with the disc scale-length (lAguerri et al. 
35 ' 



l2005UPerez et al.V2005VErwin' "20051: iMarinova & Jo gee' ^200^ 
.Laurikainen et al., 2007). Bar lengt h is also a funct ion of galaxy 
size, morphology and color ( Agu erri et al.L l2009l: iHoyle et al.L 

[mil). 

The bar st rength has b ee n det ermined by measuring 
bar torques ([ Buta & Block', 2001), isoph otal ell i pticity 



Martinet & Frie dli. 1997; Aguerri, 1999; Wh yteetaLl. 12002 : 
Marinov a & Jog eel 12007). or F ourier modes (|Qhta 



et al 



19901: Aguerri et al.[ 2Q00al: iLaurikainen et all 12005 : 
Athanassoula & Misi riotis . I2002h . This parameter depends 



on galaxy morphology. Bars in lenticular galaxies are generally 
weaker than in spirals (Das et al. 2003; Laurikainen et al. 2007; 
Barazza et al. 2008; Aguerri et a. 2009; Buta et al. 2010). 

The bar pattern speed, Qb, is the main kinematic observable 
and describes the dynamics of the bar. This tumbling pattern 
determines the position of the resonances in the disc and it is 
most usefully parametrised by a distance independent parame- 
ter = RcR/Rhar, where Rqr is the Lagrangian/corotation ra- 
dius, where the gravitational and centrifugal forces cancel out 
in the rest frame of the bar, and 7?bar is the bar semi-major axis. 
Therefore, bars that end near corotation (1< R < 1.4) are con- 
sidered fast, while shorter bars (R >l .4) are commonly called 
slow. If R < 1.0 then orbits are elongated perpendicular to 
the bar, and self consistent bars cannot exist in this regime 
(IContopoulosl D^80). The most reliable method for obtaining 
the location of corotation was that proposed by Tremaine & 
Weinberg (1984, hereafter TW method) which uses a set of sim- 



ple kinematic measurements to derive the bar pattern speed as- 
suming that the tracer obeys the continuity equation, that the 
discs are flat and that there is one well defined pattern speed. 
However, large integration times are required in medium- size 
telescopes to reach the high signal-to-noise required to apply 
the TW method. This limits its application to a small number 
of candidates. Despite the difficulties in obtaining bar pattern 
speeds, a reasonable number of nearby g alaxies have b e en in- 
vestigated (Merrifield & Kuiiken, 1995; Debattista efaP. l2002l: 
Aguerri et al., 2003; Corsini et al... .2007) finding that all bars end 
near corotation. Some of these assumptions are not applicable 
for galaxies with nested bars, and there is now a simple exten- 
sion of the TW method to multiple pattern spe eds ( Maciejewskij 
120061: ICorsini et al.l l2003l: iMeidt et all. l2009l) and the fact that 
some authors have shown that the TW metho d can be ap- 
pHed to CO (Rand & Wallin, 2004; Zimmer et 'aD.l2004l) and H^ 



velocity fie lds (Hernandez et al., 2005; Emsellem et al.L 12006 : 



iFathi et all. [2007 ; Chemin & Hernan dez, 2009: Gabbasov et al 



20091: lFathietalll2009l) opens a new window to these studies. 

Some indirect ways to derive the bar pattern speed in- 
clude methods based on numerical modelling: generating ei- 
ther self-consistent models or models using potentia ls de- 



rived from the light distribu tions ([Duv al & Atha nassoulal 



s ( [Duv 

Lindblad e t al.. 1996; Laineet^T998; WeineretaU 



1983 



2001 



Perez et al.l.l2004i: >Zanmar Sanchez et al.i. 2008) and then match- 



ing numerical experiments with the observed velocity fields; 
or by matching numerica l simulations to the galaxy ni o rphol- 
ogy dHunter et al.[. [198 8': 'Engl and [T98|jLaine et aP. [1998[: 
[Ague rri et al. [200 1[: [Rautiaine n et al.L [2005[) . Other indirect 
methods to derive the bar pattern speed include identifying 
morphological or k inemati c features with resonanc es: using a 
variety of feat ures ([Elmegre en & El megreenl [I990l) : the shape 
of dust lanes ([Athanassoula il992<) : the sign inversion of the 



1993 ): 



1995[); 



radial streaming motion acro ss co r otation (ICanziad. 
rings as resonance indicators (' ButaL [1986[: [Buta et al.[. 
phase- shift between the potential and density wave patterns 
(Zhang & Buta, 2007); lo cation of min i mum o f star formation 
(Cep a & BeckmanL i199Q[: [Aguerri et all. [2000a[) : or comparison 
of the behaviou r of the phase Fourier angle in blu e and near- 
infrared images ('Puerari & Dot toriL [1997[: [Aguerri et al.l. 119981) . 
Although possibly the most accurate indirect method to calculate 
pattern speeds is the the comparison of gas velocities to those 
obtained in numerical simulations that use a potential obtained 
from optical or near-infrared light, it is also very time consuming 
and can only be applied to a relatively small number of objects. 

The technique to determine the bar pattern speed based on 
connecting the location of rings to orbital resonances was in- 
troduced by Buta (1986). It is based on the theoretical work 
presen ted by Schwarz in a series of papers (Schwari [ T98ll 
[I984bl[ah showing how these ring structures appear near the dy- 
namical Lindblad resonances due to a bar-like perturbation. To 
directly apply this method to find the specific value of the pat- 
tern speed not only the location of the ring and the association to 
a resonance is required, but some kinematic information is also 
needed. However, we can use the ^ parametrisation of the bar 
introduced previously, and determine the ratio between the outer 
ring radius (linked to the outer Lindblad resonance, OLR) and 
the bar length. In this way, we can indirectly determine, not the 
pattern speed, but whether the bars measured are in the slow or 
fast regime. 

The bar parameters discussed above have been analysed in 
local galaxy samples. There are no previous studies in the lit- 
erature about the evolution of the length, strength and pattern 
speed of bars. In this article, we study for the first time, the dy- 
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namical evolution of bars with redshift, going from the local 
Universe to z ~ 0.8. We use a well selected sample of barred 
galaxies with outer rings to exploit the power of this method. 
The study of the dynamical evolution of bars is critical to con- 
strain the angular momentum exchange between the disc and 
the halo and their subsequent evolution. Weinberg (1985) pre- 
dicted that a bar would lose angular momentum due to a massive 
dark matter halo through dynamical friction, slowing down in 
the process. This prediction was fur t her confirmed in numerical 
simulations toebattista & SellwoodL [T998[ I200QI: IXthanassoulaL 
l2QQ3l:ISellw ood & Debattista, 200^ where they found that bars 
are slowed efficiently if a substantial density of dark matter is 
present in the region of the bar. On the other hand, if the mass 
distribution is dominated by the stellar disc, then the bar remains 
rapidly rotating for a long time. We show in this work that bars 
do not show a systematic change in their dynamical state in the 
last ~ 7 Gyrs. 

The article is organised as follows: we present the sample 
selection and morphology discussion in Sect. [J] We describe the 
method followed to measure the ring and bar radius in Sect. O 
The results are presented in Sect.|4]and we discuss their implica- 
tions in Sect. [51 Conclusions are provided in Sect.[6j Throughout 
the paper the cosmological parameters used are: //q = 70 km s"^ 
Mpc-\ Qa = 0.7, and = 0.3. 

2. Sample selection 

The galaxy samples studied in this article were extracted from 
two different surveys: low redshift galaxies were taken from the 
Sloan Digital Sky Survey (SDSS; 0.01 < z < 0.04), and high 
redshift galaxies were selected from COSMOS (0.125 < z < 
0.75). 

Two caveats must be discussed before the samples are de- 
scribed in detail: first, it is worth noticing that the galaxy samples 
are not meant to be complete in any sense, however, the selec- 
tion criteria make the two samples fully comparable. Second, in 
Aguerri et al. (2009) we studied how the resolution of the SDSS 
images can affect our detection of bars. We worked out, using 
artificial galaxies, that the shortest bars that we are able to re- 
solve have a length of ~ 9 pixels. Considering a mean PSF in 
our SDSS images with a FWHM of ir09 (2.77 pixel), we con- 
clude that we resolve bars larger than ~ 3x FWHM, or equiv- 
alently, ~ 0.5 kpc at z = 0.01 and ~ 2 kpc at z = 0.04. The 
COSMOS sample was selected using the ACS data in the F8 14W 
filter. The images were processed to a resolution of 0".' 05 pixel" | 
with an averaged PSF F WHM of a'097 (S coviUe et al.L 120071: 
iKoekemoer et al.L l2007l) . Based on the previous considerations, 
we will resolve bars larger than ~ 3xFWHM which corresponds 
to ~ 0.6 kpc at z = 0.125 and ~ 2.2 kpc at z = 0.75, matching 
perfectly the SDSS spatial resolution in the low redshift range. 

2. 1. Outer Ring Morphological Classification 

The ring morph ological classification us ed in this study is based 
on the work of iButa & Crocked (Il99ll) . They divide the outer 
rings in three main morphological classes resembling the rings 
devel oped in numerical simulations near the OLR (Schwarz, 
Il98l|) . The first class, called R'^ , is characterised by a 180° wind- 
ing of the spiral arms with respect to the ends of a bar. The sec- 
ond type is known as an ring. It is defined by a 270° winding 
of the outer arms with respect to the bar ends, so that in two 
opposing quadrants the arm pattern is doubled. The R'^ and R2 
morphologies were predicted by Schwarz (1981) as the kind of 
patterns that would be expected near the OLR in a barred galaxy. 



The third class is referred to in Buta & Crocker (1991) as the 
R1R2 morphology, where the outer arms break not from the ends 
of the bar, but from an R'^-type ring. The existence of this com- 
bined type, which may be linke d to the po pulation of both main 
families of OLR periodic orbits dSchwarzl J9811 provides some 
of the clearest evidence of the OLR in barred galaxy morphol- 
ogy. Some examples of this classification, taken from our sam- 
ple of low and high redshift galaxies, are shown in Fig. [T] Buta 
et al. (1995) derived the distribution of intrinsic axis ratios for 
the outer rings using the Catalog of Southern Ringed Galaxies. 
They found that outer rings present in barred galaxies are intrin- 
sically elliptical with an axis ratio ~ 0.82 + 0.07, and that the 
intrinsic ellipticity varies from the R'^ (~ 0.74 ± 0.08) to the R2 
(~ 0.87 ± 0.08). The intrinsic shape of the rings plays an im- 
portant role when deprojecting distances such as the bar length 
and the ring radius, thus, more intrinsically elliptical rings will 
increase the uncertainties in the measurements. We decided to 
remove from our samples the R'^ type of rings, and keep only the 
R2 types since they are intrinsically rounder. In fa ct, their intrin- 
sic shape is very similar to that of typical discs (iFasano et al.L 
1993: Rvden. .20041 

2.2. Low redshift 

The barred ringed galaxies at low redshift were obtained from 
the galaxy sample analysed in Aguerri et al. (2009). They se- 
lected a volume limited sample of galaxies from the spectro- 
scopic catalogue of the SDSS Data Release 5 (SDSS-DR5, 
Adelman-McCarthy et al. 2007). This sample covers the red- 
shift range 0.01 < z < 0.04, down to an absolute magnitude 
of Mr < -20, and with low inclination / < 60°. The full sample 
consist of 3060 galaxies with a morphological mix of 26% ellip- 
ticals, 29% lenticulars, 20% early type spirals, and 25% late type 
spirals. Galaxies were classified in barred and unbarred systems 
by searching for absolute maxima in the ellipticity radial profiles 
of their isophotes (see Aguerri et al. 2009 for details). From the 
barred sample we visually inspected the SDSS galaxy images in 
order to look for the presence of outer rings of type R2. 

We obtained a total of 1 8 barred galaxies with suitable outer 
rings features. Table [T] shows the main properties of the bars and 
rings features measured in these galaxies. 

2.3. High redshift 

As for the low redshift sample, we have chosen a number of low 
inclination galaxies from the third release of the COSMOS HST 
survey (Scoville et al., 2007). We first downloaded all the 81 
image tiles from the COSMOS/ACS fields, observed using the 
F814W(/-band) filter, from the Multimission Archive at STScI 
(MASTQ) for visual inspection. These observations cover ~ 2 
deg^ with a pixel scale (for the drizzled data) of 0^.^05 pixel"^. 

We visually scanned the COSMOS fields to look for clearly 
ringed barred galaxies. After a preliminary list was created, we 
correlated the positions with the spectroscopic redshifts from the 
Very Large Telescope (VLT) and Magellan COSMOS spectro- 
scopic surveys (zCOSMOS Survey; Lilly et al. 2007) to search 
for candidates with reliable redshifts. We searched also for pho- 
tometric redshifts for the remaining candidates. We used pho- 
tometric redshifts determined by Faure et al. (2008) using the 



^ Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope 
Science Institute. STScI is operated by the Association of Universities 
for Research in Astronomy, Inc. under NASA contract NAS 5-26555. 
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Fig. 1: Left column: schematic view of the different ring types in barred galaxies (see Sect. 2.1). Ring type examples of barred 
galaxies at high (middle column), and low (right column) redshifts. 



Le Phare photometric redshift estimation code (Ilbert et al. 
2006), details concerning the multi- wavelength photometry can 
be found in Mobasher et al. (2007). Faure et al. (2008) used 1095 
spectroscopic redshifts from the zCOSMOS Survey (Lilly et al. 
2006) to calibrate the ground-based photometric zero points. 
Using eight bands, this method achieves a photometric redshift 
accuracy of (Taz/(1+z^)=0.031. 

From this sample, only galaxies showing type R2 ring were 
included in the sample. Finally the high redshift sample consists 
of 26 galaxies. The main properties of the bars and rings fea- 
tures are shown in Table [21 The sample covers the redshift range 
0.125 <z<0. 75. Figure [2] shows the redshift distribution. 



3. Outer ring radius and bar size definition 

Our approach to quantifying the dynamical state of the bars in 
our sample of ringed galaxies is based on the measurements of 
both the bar and ring radius. We deproject bar size in the plane 
of the galaxy using /, the galaxy inclination, and 6, the position 
angle of the galaxy component (bar or ring). We assume that 
the outer ring reflects the properties of the disc, and therefore 
that the ellipticity and position angle of the ring and disc are the 



same. In Sect. 12. II we have discussed that both components are 
intrinsically similar and possible differences will only affect our 
results by introducing a large scatter. Under this hypothesis, the 
ring radius does not need to be deprojected, since it is measured 
along the major axis of the galaxy, and the galaxy inclination can 
be derived simply by / = arccos (1 - 6ring). 



3.1. Ellipse fitting 



The low redshift sample measurements were derived by using 
the ellipticity and position angle radial profiles extracted from 
the symmetrised images. This approach allows us to clean the 
images from spurious sources. It works as follows: each image 
is rotated 180° with respect to the galaxy centre. Then, we sub- 
tract the rotated frame from the original one. The residual image 
was sigma-clipped to identify all the pixels with a number of 
counts lower than Icr, where cr is the r.m.s. of the image back- 
ground. The value of the deviant pixels was set to zero. Finally, 
the cleaned image was subtracted from the original one to get the 
symmetrised image. The ellipses were then fitted to the isophotes 
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Table 1: General properties of the SDSS sample. 



Name 


Rbar-min 


Rbar-max 


PAbar 


^ring 


PA 


Rring 


z 




(kpc) 


(kpc) 


(degrees) 




(degrees) 


(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



SDSSJ104924.86-002547.5 


6 


6 


5 


1 


121.6+0.9 





402+0 


004 


142.2+0.8 


13.0 





039 


SDSSJ102543.29+393846.9 


7 


3 


5 


1 


153.4+0.2 





075+0 


005 


74.8+2.3 


13.0 





023 


SDSSJ122529.23+471623.4 


5 


7 


4 


7 


11.6+0.1 





174+0 


005 


30.1+14.3 


10.8 





025 


SDSSJ130235.73+41 1924.1 


4 


8 


3 


7 


165.4+0.7 





186+0 


005 


45.4+1.5 


9.9 





028 


SDSSJ120732.62+324846.7 


7 


1 


4 


7 


54.4+0.5 





084+0 


007 


74.1+20.0 


11.6 





026 


SDSSJ133259.13+321913.6 


5 


7 


3 


8 


177.3+0.3 





044+0 


007 


106.7+8.6 


8.7 





035 


SDSSJ012858.63-005656.3 


8 





6 


3 


6.2+0.3 





341+0 


002 


83.6+0.4 


16.4 





018 


SDSSJ083220.43+412132.0 


3 


5 


2 


4 


52.8+0.5 





048+0 


005 


111.3+10.9 


7.9 





025 


SDSSJ083630.84+040215.6 


6 


6 


4 


5 


94.2+0.2 





164+0 


009 


37.0+3.8 


12.7 





029 


SDSSJ091426.23+360644.1 


6 


3 


5 





164.4+0.0 





163+0 


031 


156.1 + 1.0 


10.1 





022 


SDSSJ123234.57+492312.2 


5 


3 


3 


4 


176.1+0.5 





074+0 


012 


114.1 + 10.7 


8.7 





040 


SDSSJ142412. 12+350846.0 


3 


6 


2 


7 


45.3+1.0 





238+0 


006 


52.9+26.7 


9.3 





029 


SDSSJ153619.30+493428.3 


5 


3 


2 


9 


17.8+2.1 





211+0 


006 


73.1+0.7 


9.7 





038 


SDSSJ160331.62+492017.3 


8 


7 


5 


2 


70.9+0.4 





255+0 


005 


32.3+7.2 


16.9 





020 


SDSSJ172721.89+593837.6 


5 


9 


4 





171.3+0.1 





202+0 


004 


148.8+1.1 


11.4 





028 


SDSSJ123313.69+121449.2 


3 


9 


3 


1 


118.5+0.9 





023+0 


004 


65.3+11.5 


6.9 





026 


SDSSJ120609. 11-025653.2 


5 


3 


3 


9 


61.5+0.7 





149+0 


009 


17.3+4.6 


11.2 





026 


SDSSJl 1 1044.88+043039.0 


6 


8 


5 





97.6+0.2 





062+0 


007 


64.6+9.2 


12.5 





029 



NOTE. Col. (1): Galaxy name from SDSS; Col. (2): bar radius calculated using the position of the minimum el- 
lipticity; Col. (3): bar radius calculated using the position of the maximum ellipticity; Col. (4): position angle of 
the bar; Col. (5): ring ellipticity; Col. (6): position angle of the ring; Col. (7): ring radius; Col. (8): spectroscopic 
redshift from SDSS 



Table 2: General properties of the COSMOS sample 



Name 


Rbar-min 


Rbar-max 


PAbar 


^ring 


PAj-ing 


Rring 


z 




(kpc) 


(kpc) 


(degrees) 




(degrees) 


(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


812947 


5.4 


4.1 


76.3+0.1 


0.068+0.010 


120.5+14.8 


9.8 


0.125 (5) 


816960 


5.5 


3.4 


174.8+1.2 


0.059+0.035 


85.2+22.6 


7.5 


0.311 {s) 


J095928.30+020109.0 


5.3 


3.8 


63.2+0.6 


0.088+0.045 


141.3+41.8 


9.1 


0.530 (/?) 


817887 


5.6 


4.6 


116.2+2.8 


0.056+0.028 


88.4+17.5 


15.1 


0.672 {s) 


823705 


3.6 


2.7 


46.4+0.7 


0.165+0.028 


74.4+17.2 


6.3 


0.491 {s) 


824759 


3.7 


2.9 


136.1+1.8 


0.161+0.023 


158.1+7.0 


7.0 


0.751 {s) 


825492 


5.1 


3.6 


108.4+2.4 


0.152+0.057 


108.1+22.2 


8.7 


0.736 {s) 


833039 


4.3 


3.3 


100.6+1.1 


0.203+0.056 


32.8+ 4.8 


8.8 


0.360 {s) 


J100233.98+022524.3 


4.7 


2.2 


172.7+0.9 


0.167+0.015 


63.1+4.4 


9.4 


0.720 {p) 


J095938.81+020658.7 


10.3 


7.6 


111.7+0.4 


0.071+0.035 


149.6+75.9 


16.6 


0.409 {p) 


J095935.08+020127.2 


4.3 


3.0 


176.8+79.2 


0.158+0.044 


26.3+36.9 


10.3 


0.357 (/7) 


J100204.95+022739.7 


5.8 


3.7 


124.4+1.0 


0.116+0.019 


122.7+ 6.6 


10.5 


0.507 {p) 


841055 


5.7 


3.4 


134.9+0.6 


0.188+0.020 


178.0+56.2 


13.7 


0.376 {s) 


J095759.45+022810.5 


4.4 


2.6 


20.6+0.9 


0.042+0.013 


131.0+23.0 


9.2 


0.119 {s) 


851598 


6.4 


4.4 


9.6+0.5 


0.080+0.047 


102.2+30.7 


10.0 


0.346 {s) 


852495 


6.8 


5.7 


96.7 +3.1 


0.265+0.015 


21.9+ 1.5 


12.2 


0.705 {s) 


852636 


7.1 


4.9 


0.7+80.0 


0.040+0.026 


125.6+50.7 


12.5 


0.345 {s) 


852155 


7.7 


5.6 


173.0+0.3 


0.306+0.009 


160.4+ 0.8 


15.3 


0.305 {s) 


J100254.88+024645.8 


4.4 


3.2 


18.7+1.5 


0.065+0.036 


128.0+45.7 


8.5 


0.468 {p) 


840577 


5.1 


3.5 


68.3+2.4 


0.085+0.023 


89.5+15.0 


7.9 


0.539 {s) 


838743 


4.4 


2.8 


48.7+0.5 


0.144+0.025 


60.0+3.8 


8.9 


0.126 (5) 


830974 


9.6 


5.3 


11.5+1.4 


0.149+0.025 


57.7+7.0 


13.9 


0.695 {s) 


811921 


4.9 


3.1 


25.2+1.6 


0.068+0.027 


89.5+30.7 


8.2 


0.371 {s) 


813153 


7.5 


6.0 


60.0+0.7 


0.159+0.069 


47.3+28.1 


13.2 


0.529 (5) 


831775 


3.6 


1.8 


125.1+0.1 


0.140+0.028 


156.8+76.5 


6.2 


0.381 {s) 


J100217. 12+023024.1 


4.9 


3.4 


157.8+0.8 


0.083+0.022 


97.7+15.6 


8.1 


0.379 {p) 



NOTE. Col. (1): Galaxy name from COSMOS; Col. (2): bar radius calculated using the position of the minimum 
ellipticity; Col. (3): bar radius calculated using the position of the maximum ellipticity; Col. (4): position angle 
of the bar; Col. (5): ring ellipticity; Col. (6): position angle of the ring; Col. (7): ring radius; Col. (8): {s): 
spectroscopic redshift from zCOSMOS; (/?): photometric redshift from Faure et al. (2008). 



of the symmetrised images of the galaxies using the IRAlQ task ELLIPSE ("Jedrzejewski, "1987). We used an iterative wrapped 

~; \ ~ :; ~~ , . , . , , procedure which runs the ellipse fitting several times, changing 

2 IRAF is distributed by NOAO, which is operated by AURA Inc., ^ f & ' & & 

under contract with the National Science Foundation. 
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Fig. 2: Redshift distribution for the high redshift sample. 



the trial values at each fit iteration, until a good fit at all radii 
is obtained. At each fixed semi-major axis length, the coordi- 
nates of the centre of the fitting ellipse were kept fixed. This 
centre was identified with the position of the central intensity 
peak. The trial values for the ellipticity and position angle were 
randomly chosen between and 1 and between -90° and 90°, 
respectively. The fitting procedure stopped when either conver- 
gence was reached or after 100 iterations. 

The high redshift sample ellipticity and position angle pro- 
files were derived using the same wrapping procedure to max- 
imise the goodness of the ellipse fitting. However, in this case we 
preferred not to symmetrise the images but apply a 2 x 2 pixels 
box smoothing. The image symmetrisation was not needed since 
the sample galaxies were not contaminated by other sources 
within their projected surface and the smoothing provided bet- 
ter radial profiles by improving the signal-to-noise ratio of the 
images. 

3.2. Bar length 

The ellipticity radial profile of a typical barred galaxy reaches 
a minimum at the centre, because of either seeing effects or a 
spherical bulge. Then, it usually shows a global increase to a lo- 
cal maximum, and then suddenly decreases towards a minimum 
at the location where the isophotes become axisymmetric in the 
face-on case. The position angle profile is also characteristic in 
barred galaxies, being constant in the bar r egion and then chang- 
ing to fit the outer d isc orientation (e.g.. IWozniak et all 1 19951: 
I Aguerri et al.Ll200()bl) . These characteristic profiles are produced 
by the shape and orientation of t he stellar orb its of the bar (see 
IContopoulos & GrosbolL 119891: lAthanassoulaL 119921) . Different 
methods have been used to measure the bar lengt h based on the 
ellipt i city and position angle radia l profi les (see lAthanassoulal 
I20Q2I: iMiche l-Dansac & Wozni"^ 12006). However, the solu- 
tion is always ambiguous and it can lead to misleading re- 
sults. In order to remove these uncertainties we decided to 
measure the bar length as the midpoint between the radius 
of the maximum and minimum ellipticity. These two differ- 
ent measurements of the bar length represent the extreme cases 
iMichel-Dansac & Wozniak, 2006) and therefore they represent 
an upper limit of our errors in the bar radius measurements. We 
preferred this solution not to bias our conclusions. The position 



angle of the bar, which is needed to deproject the bar length, was 
measured at the position of the maximum ellipticity, so we avoid 
problems related to position angle variation in the bar-disc re- 
gion. An example of this method applied to three of our low and 
high redshift galaxies is shown in Figs.[3l|4l respectively. 

3.3. Ring radius 

As for the bar component, the ring radius was derived based on 
the ellipticity and position angle radial profiles. In the ring re- 
gion, we expect that the ellipticity and position angle radial pro- 
files will remain constant due to the stellar orbits in the ring. 
Therefore, we identify the region of the profile where the ring is 
present and we measured the ring radius as the position where 
the ellipticity and position angle become constants. The ring el- 
lipticity and position angle also needed for deprojecting, were 
derived as a mean of these constant values. The error in the ring 
radius has been calculated by comparing the estimated ring ra- 
dius with the radius at which the ellipticity varies more than 3 
times the standard deviation of the disc ellipticity. Figures lA.ll 
and I A. 21 from the Appendix show all the galaxies with the ring 
radius overplotted. 



4. Results 

Table [T] and Table [2] shows the obtained parameters for the ring 
radius, ellipticity, position angle and the bar semi-major axis, as 
derived in Sect. 13. 11 Most galaxy inclinations lie below / < 40°. 
The bar size range, using the maximum ellipticity, covers from 
2.5 to 6.3 kpc. Most of the bars in the local Universe (about 70%, 
see Aguerri et al. 2009) are within this bar size range. Similar 
values of the bar size range are found for our high redshift galax- 
ies. The mean bar radius of our low and high redshift galax- 
ies are 4.5+1.04 and 3.5+1.33 kpc, respectively. This means 
that within the errors both galaxy samples have similar bars ac- 
cording to their lengths and similar to local samples of barred 
galaxies (see Aguerri et al. 2009). The average bar size, using 
the minimum ellipticity, for our low and high redshift galaxies 
are also similar: 5.9+1.43 and 5.3+1.70 kpc, respectively. Thus, 
both samples of galaxies show similar bar sizes independent 
of the m ethod used for determining the b ar length. It has been 
argued (Michel-Dansac & Wozniak, 2006) that the sizes calcu- 
lated using the minimum of the ellipticity correlate well with the 
position of corotation, giving a more physically significant size 
than measurements obtained with the maximum of the ellipticity, 
which clearly underestimates the true bar size. To avoid prob- 
lems related to the bar size calculation, as explained in Sect. 13.21 
we have opted for using the mid-point and to take into account 
the values of Rbar using both methods to obtain the errors. 

We have determined the strength of the bars for the low and 
high redshift galaxies by using the maximum ellipticity of the 
bar (see Aguerri et al. 2009). Both samples cover the same range 
of bar strengths. Thus, the mean values of the bar strength of our 
low and high redshift samples are: 0.20+0.07, and 0.17+0.05. 
These values are similar to the mean strength of bars in the local 
Universe (0.20+0.07; see Aguerri et al. 2009). We can conclude 
that according to the size and strength of the bars, our low and 
high redshift galaxy samples have similar bars as those found in 
a complete local sample of barred galaxies (see Aguerri et al. 
2009). 

To determine whether our galaxies are in the fast or slow 
range (see Sect[T]) we define the ratio 'Rnng= Rring/Rbar, where 
Rring is the ring radius and R^^^, is the bar semi-major axis, as 
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characterised in Sect. [3] Because we cover this ratio for galax- 
ies with redshifts between < z < 0.8, we can study possible 
changes of this ratio with redshift. Figure [5] shows the distribu- 
tion of ring radii (Rring) vs. the bar semi-major axis (Rbar) for 
the whole sample. We consider a fast bar those bars for which 
the RcR/Rbar T^tio lics between 1.0 and 1.4. This ratio has been 
plotted in Fig. [5] for both values and is calculated using linear 
resonance theory and a flat rotation curve (Athanassoula et al. 
1982), in this case the position of the OLR (i.e., the ring radius) 
and the CR are related in the following way: 



(1) 



where 6 lies between 0.7 and 1.0 for early type discs, see 
Athanassoula et al. (1982). We take fcl.O in Fig.Ofor simplicity 
but this choice does not alter the results. It is clear from Fig. [5] 
that all the galaxies, independent of their redshift bin, fall into 
the fast-bar category. 

We have investigated the influence of the inclination of the 
galaxies on this result. Thus, the average values for at dif- 
ferent inclinations are the following: i) for the low redshift sam- 
ple and b/a > 0.9, the average ^nng = 0.51+ 0.06, for b/a < 0.9, 
the average %mg = 0.50 ± 0.08 ii) for the high-redshift sample; 
for b/a > 0.9, the average ^ring = 0.55 ± 0.09 and for b/a < 0.9, 
the average "Rnng = 0.52 ±0.12. All the values, independently of 
redshift and inclination are comparable; and therefore, we do not 
see changes of this ratio with redshift. 



Our low and high redshift samples are by no means com- 
plete. Therefore, it is customary to investigate whether this fact 
is affecting the results presented in Fig. [5J Since our low and 
high redshift barred galaxies are similar in size and strength we 
could be biasing the resulting pattern speeds toward a particular 
regime. In other words, the fact that we have not observed evo- 
lution in the pattern speed could be just due to the selection of 
similar fast bars. However, we know from a study of a complete 
sample of local barred galaxies (Aguerri et al. 2009), that only 
30% of the local bars show larger lengths than our r inged barred 
galax i es. Studies of high re dshift bars, 0.4 < z >0.8, (I Jo gee et al.L 
l2004l: iBarazza et al.L l2009l) have shown that the bar size distri- 
bution is similar to that of local galaxies and; therefore, as dis- 
cussed before, similar to the bars size range of the bars presented 
in this work. Therefore, we do not seem to be looking at any 
special type of bar by analysing ringed galaxies. As previously 
mentioned, from numerical models, bars get longer and slower 
as they age. We can then set a 30% upper limit to the bars that 
could have suffered a change in their pattern speed in the last 
7 Gyr, assuming that nearby long bars are the end-products of 
the evolution of fast bars. The remaining 70% of bars did not 
substantially lose angular momentum to the halo, maintaining 
their pattern speed. This discussion might be related to an im- 
plicit morphological bias, since it remains, even for local galax- 
ies, to derive the pattern speed of bars in very late-type gaseous 
rich spirals which might suffer an intrinsically different evolu- 
tion (e.g., Bournaud and Combes 2002). 



5. Discussion 

5. 1. Possible caveats 

The tightness of the results shown in Fig. \5\ is somewhat un- 
expected considering the intrinsic uncertainties inherent to the 
measurements used in this work. For instance, we have assumed 
that the outer ring are perfectly circular, which is critical for the 
deprojection of the ring and bar lengths. In Sect. [3] we justified 
the assumption of roundness for the R2 and R1R2 ring morpholo- 
gies. In addition, we have taken int o acc ount the limits of intrin- 
sic axis ratios given m iButa et al.1 (fl^5 ) to calculate the errors 
in the projected sizes, and we have demonstrated in the previous 
section that our result does not depend on the inclination of the 
galaxies. Therefore we conclude that projection effects are not 
biasing our results. 

Another possible caveat to our result might be that the 
choice of ring galaxies biases the sample towards a certain pat- 
tern speed domain. However, numerical simulations (iByrd et al.L 
[T994) have shown that resonant outer rings can be present in 
both fast and slow bars. All types of ring morphologies are 
found at different pattern speeds. In the same work, all types 
of ring morphologies also developed for different bar strengths. 
Athanassoula et al. (2010), who presented a new theory for ring 
and spiral formation, argues that there is a connection between 
the bar strength and the morphology of the rings. Nevertheless, 
R2 rings, as are those selected in this work, are located in 
barred galaxies with similar bar strengths as our galaxies (see 
Athanassoula et al. 2010). In addition, they show that R2 type 
rings can be formed in galaxies with fast and slow bars. 

Although the bar size as measured in different rest-frame 
band passes could be different, it has been recently shown 
(Gadotti 201 1) that in fact the difference in bar size is negligible 
and therefore we are not introducing a bias by measuring the bars 
a high-redshift near the g-band rest- frame while the low-redshift 
bar sizes are measured from the SDSS r-band. 



5.2. Comparison with the results from numerical 
modelling 

A recent numerical work (I Villa- Vargas et al.1. l2010l) shows that 
the evolution of the pattern speed and bar-growth of a bar em- 
bedded in a live dark matter halo depends strongly on the gas 
content. In their simulations, a fixed fraction of the total mass 
was converted to gas mass, and the evolution of the bar param- 
eters is then followed in time. The presence of gas changes the 
evolution of both the bar growth and the pattern speed evolution, 
the addition of gas can stop, or even speed-up, the pattern speed 
of the bar with time. The bar size is anti-correlated with the disc 
gas fractions. These gas-rich galaxies would be related to early- 
type galaxies because the gas leads to larger central mass con- 
centration and therefore larger bulges. The results we present 
in this paper could be in agreement with these gas-rich models. 
However, the full picture is still unclear since it is observed that 
longer bars reside in late-type galaxies (e.g., Erwin 2005) which 
is against the model predictions. Furthermore, we should then 
explain why all the galaxies should have started with similar gas 
fractions in their discs. 

There is also the possibility that the bars that we see at 
z ~ 0.8 do not survive till the present and therefore, we do 
not see evolution because the time- scales involved in the for- 
matio n and destruction are too s hort. It has been discussed 
(see iPfenniger & N ormanL 1 19901: iBournaud & CombesL l2002l: 
iBournaud et al.L 12005) that gas-rich bars, i.e., late-type spirals, 
are short lived, with lifetimes of 1-2 Gyr. This short time scale 
would mimic a lack of evolution of long-lived bars; however, 
the galaxies in the sample show morphologies typical of early 
type spirals and there is evidence, from stellar popu lation studies 
dPerez et all. l2009l: ISanchez-Blazquez et al.1 lioilh . that bars in 
early type galaxies are long-lived. If this is the case, and most of 
the ring galaxies we observe present long-lasting bars, it would 
imply that bars cannot have grown in time and kept being in 
the fast speed regime without increasing significantly in size. 
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arcsec R (arcsec) R (arcsec) arcsec 

Fig. 3: Three examples of the ringed galaxies selected from the SDSS data. From left to right: original r-band SDSS image, ellipticity 
isophotal radial profile, position angle radial profile of the isophotes, and deprojected image. The vertical full, dashed and dotted- 
dashed lines represent the Rbar-min, Rbar-max, and Rring, respectively. The horizontal full, dashed and dotted-dashed lines shows the 
PA of the bar, ellipticity of the disc and PA of the disc, respectively. The circle represented in the right-most panels has a radius 
equal to the measured radius of the ring. 



Therefore, the fact that we see the ring radius and the bar size 
covering the same size range at all redshift, and moreover large 
bars at high-redshift, implies that bars do not grow significantly 
in size with time. 



The result shown in Fig. [5] implies that bars have not evolved 
considerably, neither in size nor in pattern speed, since around 
the time when the Universe was half its present age. Most numer- 
ical simulations obtain bars that evolve with time, getting longer 
and stronger while slowing down (e.g., Debattista &Sellwood 
1998; Athanassoula 2003). This effect is mostly due to the an- 
gular momentum exchange of the bar-disc system with the dark 
matter halo. Thus, the fact that bars are compatible with fast ro- 
tators at all redshifts indicates that the angular momentum ex- 
change between the bar and halo has not been important enough 
in the last 7 Gyr to slow down bars. If the pattern speed can be 
used to set constrains to the halo-to-disc mass ratio, these results 
might imply that the discs in the high surface-brightness galaxies 
of our sample are maximal. 



6. Summary and conclusions 

We have analysed 44 low inclination ringed galaxies spanning a 
redshift range between < z < 0.8 to study the possible evolu- 
tion of the pattern speed in the last 7 Gyrs. 

We calculated for each galaxy a morphological parameter 
indicative of the dynamical state of their bars. In particular, we 
derived whether they are fast or slow rotators. We find that the 
bar pattern speed does not seem to change with redshift and that 
all bars are compatible with being fast bars. 

If the bars analysed are long-lasting, their size and bar 
strength have not significantly changed in time. The fact that, 
independent of the redshift, the bars are fast rotators and their 
size has not significantly changed in time could have also large 
implication for bar evolution models that mostly predict a bar 
growth with time. It has been argued that the exchange of an- 
gular momentum with a centrally dense halo causes the bar to 
evolve; however the present results might imply that the disc 
in the high surface-brightness galaxies is maximal and the cen- 
tral mass density is dominated by the stellar component which 
would lower the angular momentum exchange between the disc 
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Fig. 4: Three examples of the ringed galaxies selected from the COSMOS data. From left to right: original r-band COSMOS image, 
ellipticity isophotal radial profile, position angle radial profile of the isophotes, and deprojected image. The vertical full, dashed 
and dotted-dashed lines represent the Rbar-min, Rbar-max, and Rring, respectively. The horizontal full, dashed and dotted-dashed lines 
show^s the PA of the bar, ellipticity of the disc and PA of the disc, respectively. The circle represented in the right-most panels has a 
radius equal to the measured radius of the ring. 



and the halo and slow down the bar evolution (e.g., Debattista & 
Sellwood, 2000, but for a different conclusion, see Athanassoula, 
2003). 

This is the first time that the pattern speed evolution has been 
investigated from the observational point of view. The results 
presented here place strong constrains on the bar evolution mod- 
els. 
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Fig. A.l: Low redshift sample with black solid circles indicating the ring size. 
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Fig. A.l: Low redshift sample with black solid circles indicating the ring size. 
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Fig. A.2: High redshift sample with black solid circles indicating the ring size. 
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Fig. A.2: High redshift sample with black solid circles indicating the ring size. 
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Fig. A.2: High redshift sample with black solid circles indicating the ring size. 
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